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Abstract 

Calculation of the vacuum polarization, ((f> 2 (x)}, and expectation value of the stress tensor, (T M „(:r)), has seen a recent 
resurgence, notably for black hole spacetimes. To date, most calculations of this type have been done only in four dimensions. 
Extending these calculations to d dimensions includes d-dimensional renormalization. Typically, the renormalizing terms are 
found from Christensen's covariant point splitting method for the DeWitt-Schwinger expansion. However, some manipulation 
is required to put the correct terms into a form that is compatible with problems of the vacuum polarization type. Here, after 
a review of the current state of affairs for {<f> 2 (x)} and {T^{x)) calculations and a thorough introduction to the method of 
calculating (4> 2 (x)}, a compact expression for the DeWitt-Schwinger renormalization terms suitable for use in even-dimensional 
spacetimes is derived. This formula should be useful for calculations of {(f> 2 (x)) and (T M!/ (a;)) in even dimensions, and the 
renormalization terms are shown explicitly for four and six dimensions. Furthermore, use of the finite terms of the DeWitt- 
Schwinger expansion as an approximation to (4> 2 (x)} for certain spacetimes is discussed, with application to four and five 
dimensions. 
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I. INTRODUCTION 

A. Vacua and particle creation 

Hawking radiation shows that a black hole and its horizon are spacetime regions where gravitational effects become 
important in quantum field theory. The initial approach followed by Hawking [l| was to study quantized ingoing and 
outgoing field modes in a fixed black hole background. This process yields a flux of particles, produced and emitted 
out of the surrounding quantum vacuum with a thermal spectrum. In a flat spacetime the vacuum is uniquely defined 
and relatively easy to find, but in a curved spacetime it is not so simple. There are cases where it may be impossible 
to define a vacuum state, as in strongly time dependent geometries. On the contrary, there are situations where many 
vacua exist for a given geometry. The Schwarzschild solution, for example, has three possible vacua: Boulware 0], 
Unruh [3], and Hartle-Hawking jj], all three of which are interesting and important. The Boulware vacuum describes 
the vacuum in a region near the surface of a highly compact star, where it has a small negative energy density. Near 
the horizon of a black hole, however, this energy density blows up, so the Boulware vacuum is inappropriate for a black 
hole geometry. On the other hand, the Unruh and Hartle-Hawking vacua are both consistent for the Schwarzschild 
black hole. The Hartle-Hawking vacuum has a small, finite, negative energy at the event horizon, which in turn is 
responsible for the production of particles and subsequent Hawking radiation at a given temperature. The Unruh 
vacuum, also having a small, finite, negative energy at the event horizon, is produced by the complete gravitational 
collapse of an object, and in this sense is more physical. However, the one that usually simplifies the calculations is 
the Hartle-Hawking vacuum, and, if desired, it is possible to pass from this to the Unruh vacuum via appropriate 
transformations. 

Once a consistent vacuum for the geometry in question, such as a black hole spacetime, has been defined, the 
associated quantities of interest may be found. For a given field <f>{x) and vacuum |0), where x represents a spacetime 
point, these quantities may be the vacuum expectation value of the field operator 4> 2 {x) i.e. (0|</> 2 (:c)|0), or (<fr 2 (x)} for 
short, and its associated vacuum expectation value of the stress-energy tensor T^ v (x) i.e. (0|T (ttI/ (a;)|0), or {T^ix)) 
for short. The quantity (cf) 2 (x)) is a useful tool in the study of quantum effects in curved spacetimes. When properly 
renormalized it gives information about vacuum polarization effects and spontaneous symmetry breaking phenomena, 
although since it is a scalar it does not distinguish between future and past surfaces, such as horizons and infinities. 
The quantity {T^^x)) provides information about the energy density and particle production. Moreover, since in 
general relativity Einstein's equations relate the spacetime curvature to the distribution of matter as encoded in 
the stress-energy tensor, for quantum fields the expectation value {T^ v (x)) is used to determine how the underlying 
geometry responds to suitable averages of the quantum fields. These back-reaction effects by the quantum fields on 
the background spacetime are described by the semiclassical Einstein equations — 8Tr(T flI/ (x)) (we use G = 1, 
c=l, n=l) (see (ESQ 

for reviews and careful explanations). The problem of quantum back reaction is certainly 
significant in the case of black holes and in other spacetimes with horizons, such as de Sitter spacetime. For instance, 
in the black hole case it leads to the complete evaporation of the black hole. 

B. Renormalizing the vacuum 

Since (4> 2 (x)) and {T^ v {x)) are constructed with product expressions, bilinear in the field operators and evaluated 
at the same spacetime point, the vacuum expectation value of these quantities diverges. For a theory to have any 
physical meaning it must give finite results, so some process must be employed to render these quantities finite. Such 
a process amounts to subtracting off some "unphysical" infinite terms. In flat spacetime, standard normal ordering 
techniques and other procedures in quantum field theory work well in regularizing and renormalizing the fields, where 
regularization identifies the infinities, and renormalization eliminates them. In general relativity, however, the energy 
density itself is a source of curvature. Therefore, when working with a quantum field theory whose energy density 
is formally divergent we must be very careful about what may be dismissed as unphysical. The standard techniques 
used in flat spacetimes do not work in curved spacetimes, so one of the great difficulties in understanding quantum 
processes in a black hole - or any other curved - spacetime is the implementation of consistent regularization and 
renormalization schemes. 

Fortunately, there are severalgenerally accepted consistent, covariant regularization and renormalization schemes 
for curved spacetimes (see e.g. Q). Of these, the most widely and consistently used method for (</> 2 (x)) and (T Ml/ (a;)) 
calculations is that of isolating the divergent terms of the DeWitt-Schwinger expansion. This covariant geodesic point 
separation method, developed by Schwinger [1], DeWitt 0, [l(|, and Christensen [lj, [l^, EH (see Barvinsky and 
Vilkovisky in Ref. [lj| for further information), is now usually called the point splitting method. The idea of the point 
splitting method is that the operator in each product is moved along a geodesic to a nearby spacetime point. The 
point separated object is expressed in terms of Green's functions G(x,x'). In the coincidence limit, where the nearby 
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spacetime point x' approaches the original point x, there will be terms diverging logarithmically (in even dimensions) 
and as inverse powers of the point separation. This point splitting method leads naturally to the DeWitt-Schwinger 
expansion, which gives an approximation for the Green's function G(x, x') when the points x and x' are separated by 
a small geodesic distance, s, along the shortest geodesic connecting them. The result is actually expanded in powers 
of the field mass m, with expansion coefficients a.k expressed in terms of geometrical quantities constructed from the 
Riemann tensor. Other renormalization methods exist, for example, dimensional continuation flEI ]. Many of these 
methods have been shown to be equivalent to the DeWitt-Schwinger approach [Icj |. 

The divergent terms of the DeWitt-Schwinger expansion are then the renormalizing counter terms to be subtracted 
from the unrenormalized exact expression for G(x,x'), prior to taking the limit x — > x' . Now the expressions 
for ((f> 2 (x)} and (T m „(ie)), written conveniently in terms of Green's functions, are in fact the properly renormalized 
expressions ((f> 2 (x)} rcn and (T^(x)) lcn . These renormalized values are those which provide information on spontaneous 
symmetry breaking and particle production, as well as being essential for calculating the backreaction by quantum 
fields on the spacetime. This means that semiclassical general relativity has physical meaning when described by 
the equations G M „ = 8Tr(T^(x)) rcn . With this renormalization process in hand, a complete set of mode functions 
and their associated creation and annihilation operators must be found by solving the field equations for (j)(x). The 
vacuum expectation values (<^> 2 (a;)) re n and (T MJ/ (a:)) ren are then put in terms of the field operators. In general the 
result is a sum over products of mode functions and their derivatives. The sum can, at least in principle, be performed 
and a finite result is achieved. 



C. DeWitt-Schwinger estimates for {(f> 2 (x)) 



There is an additional pay-off when using the DeWitt-Schwinger expansion with Christensen's point separation 
method. Since the expansion is in inverse powers of the field mass m, it is valid for many spacetimes provided m 
is large enough. In this case the finite terms of the expansion can provide approximations for both ((f> 2 (x)) len and 
(T^v{x)) Icn . For instance, given a scalar field 4>{x) the Feynman Green's function corresponds to {(jP{x)) (see e.g. 
Q), so that the finite terms of the expansion directly give a physical (0 2 (x)) ren . The major obstacle in the DeWitt- 
Schwinger expansion is to compute the coefficients afc. For a scalar field the first three coefficients, ao, a\, and a2, 
have been computed by DeWitt [§, [l(| ; the coefficient 03 has been computed in the coincidence limit by Gilkey ; 
and, the coefficient 04 has been computed in the coincidence limit by Avramidi [17| . and by Amsterdamski, Berkin, 
and, O'Connor [HI]. Additionally, Barvinsky et. al. [l^| have calculated these coefficients using different methods. 
Thus, Christensen's method plays definitely two roles - it is the basis for point splitting renormalization, and it yields 
an estimate for the quantities (4> 2 (x)) Yen and (T^ v (x)) rcn . 

It should be stressed that the DeWitt-Schwinger expansion together with Christensen's point separation method is 
an approximation that does not hold in all regions of all spacetimes. For example, the results of Kay and Wald (20j 
show that for a Reissner-Nordstrom black hole in asymptotically de Sitter spacetime, (<fi 2 (x)) cannot be regular on 
both the event and cosmological horizons when these horizons have unequal temperatures. 



D. Calculations of (4> 2 (x)) and applications 

Christensen's work is quite general, and in principle can be applied to any spacetime. For cosmological as well 
as some black hole applications, see Ref. [5j for works up to around 1980. Many other examples can be given. 
For massless scalar fields analytical results were reported by several authors. Candelas studied a massless scalar 
field minimally coupled in the Schwarzschild geometry, where (4> 2 (x)) and (T^ v {x)) were worked out on the event 
horizon [2l|. Candelas and Howard [22j and Fawcett and Whiting [23j] extended the calculation of (cf> 2 (x)) to the 
exterior region. Candelas and Jensen [24| calculated {4> 2 (x)) in the interior region, and finally Howard and Candelas 
[25L l26j and Fawcett [13] calculated {T^ v {x)) for the whole of Schwarzschild, definitively extending the pioneering 
work of Candelas [2l[. In this context of Schwarzschild black holes it was shown by Hawking [28| and Fawcett 
and Whiting (23[ that the mean square field ((f> 2 (x)} can give considerable insight into the physical content of the 
different possible vacua and in the study of theories with spontaneous symmetry breaking. Massless scalar fields 
in Reissner-Nordstrom and Kerr-Newman spacetimes were studied by Frolov |29j, where (4> 2 (x)} was found on the 
event horizon of a Reissner-Nordstrom black hole and on the pole of the event horizon of a Kerr-Newman black hole. 
Since analytical and numerical calculations are difficult, approximation schemes have been devised for calculating 
{4> 2 {x)) and (Tfj,v(x)) for Schwarzschild, Reissner-Nordstrom, and Kerr-Newman black holes [3(1 HH, HH, HI; H3 ■ For 
massless electromagnetic fields several works have calculated (T^(x)) for Schwarzschild, Reissner-Nordstrom, and 
Kerr-Newman black holes [35|,[36|, H3, 38]. All these works are for massless fields, where calculations may simplify due 
to the conformal invariance of the system. 
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A new approach came with the work of Anderson [39|, Ecj where the method was applied consistently to massive 
scalar fields. In Ref. [3§| (cj) 2 (x)) was calculated for a generically coupled massive scalar field in the Schwarzschild 
geometry, and in [4(3] a powerful formalism was laid down for finding (cj) 2 (x)) in a general static spherical geometry, 
which includes Schwarzschild and Reissner-Nordstrom solutions. This was possible through the use of the Plana sum 
formula which converts sums into integrals. The method was extended by Anderson, Hiscock, and Samuel 4 lL_ 4211 to 
find (T^ u (x)) rcn for a massive scalar field in a general, static, spherical geometry. The approach of Refs. [40l l4ll. |4^| 
uses a Wentzel-Kramers-Brillouin (WKB) approximation for the mode functions to compute (<j) 2 (x)) and {T^ v (x)) 
to orders rn~ 4 and m~ 2 respectively. It is further found that, when applied to the Reissner-Nordstrom spacetime, 
the DeWitt-Schwinger expansion provides values quite close to the numerical results when the field mass mM > 1, 
where M is the black hole mass (we put h = 1). Anderson's approach has been developed and applied to other 
cases. Cylindrical black hole spacetimes have been examined by DeBenedicti s 1431 . who worked out {<p 2 {x)) for scalar 
fields, and by Piedra and Oca, who have studied spinor fields [44[. Sushkov j45| studied it for wormholes, Berej and 
Matyjasek [46| for the spacetime of a nonlinear black hole, Satz, Mazzitelli, and Alvarez 1471 f or the vacuum outside 
stars, and Winstanley and Young [48| for lukewarm black holes. Finally, Flachi and Tanaka[49l] have used Anderson's 
method to compute (<f> 2 (x)} in asymptotically anti-de Sitter black hole geometries. New approaches have been devised 
by Anderson, Mottola, and Vaulin [50(, whereas Popov and Zaslavskii have discussed the WKB approximation in the 
massless limit 1511. 



E. Renormalization in d dimensions and this paper 

Christensen has remarked [HI that while his methods are valid in arbitrary dimensions, where the procedure is 
the same as in four dimensions, calculating quantities such as {(p 2 (x)) and {T^{x)) in higher dimensions "would be 
extremely long and would probably have to be done on a computer," mainly due to the complexity of the renor- 
malization problem. This comment is still true and the very few works since 1978 that have tried to come to terms 
with the renormalization techniques in curved <i-dimensional spacetimes do prove the difficulty of the extension of 
the procedure. Nevertheless the interest in these techniques to spacetimes with more than four dimensions has been 
renewed as the result of progress in areas such as string theory, AdS/CFT (anti-de Sitter/conformal field theory) 
conjecture, Kaluza-Klein theories, extra large-dimensional scenarios, and the related brane world scenarios. 

Earlier, Frolov, Mazzitelli, and Paz [52| studied polarization effects in black hole spacetimes in higher dimensions. 
In the context of black holes in a braneworld, Casadio [53| discussed back reaction issues. In a very thorough work 
Decanini and Folacci [54| expressed the DeWitt-Schwinger representation of the Feynman propagator as a Hadamard 
expansion for even and odd dimensions which clearly exhibit the divergent and the regular parts of the DeWitt- 
Schwinger representation. In [55, 56] these authors presented the first explicit calculations of the stress-energy tensor 
in an arbitrary spacetime of d dimensions and provided an expression for d = 6 in the large mass limit. Following the 
ideas developed in Christensen and Fulling [53], Morgan, Thorn, Winstanley, and Young [58[ have worked out some 
properties of (T„„(a:)) for <i-dimensional spherical black holes. Herdeiro, Ribeiro, and Sampaio 59] studied the scalar 
Casimir effect on a <i-dimensional Einstein static universe where renormalization techniques are also used and where, 
incidentally, the Plana sum formula (also called the Abel-Plana formula) has been applied - in fact the formula was 
used for the first time in the context of renormalization techniques for the Casimir effect by Mamaev, Mostepanenko, 
and Starobinsky (60]; see the review [6X1 ] . 

In this paper we use the techniques developed by DeWitt [l(|, Christensen [O, [H, [II] , and Anderson [|(| EH E3] (see 
also [Il|,[22j]) and apply them to the problem of renormalization of the divergent quantity (<p 2 (x)} for a massive scalar 
field 4>(x) in a d-dimensional static spacetime, carryin g ou t the renormalization by the point splitting technique. The 
derivations presented by DeWitt [l(| and Christensen [111 [l2l [l3j are quite mathematical in nature, and the end result 
is not in a form that is amenable for renormalizing ((j) 2 (x)}. The purpose of this paper is to present a compact formula 
for the renormalization terms that may be applied to ((f> 2 (x)) calculations, which we achieve for even dimensions. As 
applications of our results in even dimensions, we single out d = 4 and d = 6. In d = 4 we compare our results 
with previous results, and surely, it is the most important dimension. We then have chosen d — 6 both because it is 
the simplest case after d — 4 and can be consistently realized if one advocates extra large dimension or braneworld 
scenarios. Odd dimensions may require other methods to find a compact formula for the renormalization terms. In 
the calculation we also find (^> 2 (a;)) r en in first approximation for large enough field masses, in both even and odd 
dimensions. We give as examples the cases d = 4 and d = 5. Again, d = 4 is singled out because it is the most 
important dimension and it can be compared immediately with the previous results of other authors, and d = 5 is 
the first odd higher dimension, and could as well be important in scenarios with large extra dimensions In brief, there 
are two purposes: one is to kill the divergences in (<j) 2 (x)} lcn ; the other is to extract the finite part of (</> 2 (x)) rcn . 

The paper is organized as follows. In Sec.[TT]the calculation of (4> 2 {x)) is thoroughly reviewed, including a discussion 
of the connection between Green's functions and operator theory, and an outline of the standard method for computing 
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(4> 2 (x)) in a static spacetime. This motivates the need for finding a compatible expression for the renormalization 
terms and shows what form they must take. In Sec. IIII Al the DeWitt-Schwinger expansion for d dimensions is 
presented, and isolation of the divergent terms is reviewed. In Sec. IIII Bl even dimensions are studied. Specifically, in 
Sec. IIII B 11 an integral representation for the modified Bessel function K v {z) in the limit of vanishing z is derived for 
even-dimensional spacetimes. For scalar fields of zero temperature, this integral representation may be used in the 
expression for the divergent terms. For a scalar field at temperature T, further manipulation is required to make the 
expression for the divergent terms useful for ((p 2 (x)) calculations. The Plana sum formula is employed to convert the 
integral into a sum plus residual terms, leading to a suitable formula for the renormalization terms in the nonzero 
temperature case. As an example, the renormalization terms are found for four- and six-dimensional spacetimes in 
Sec. IIII B~2l Section ITVl discusses estimating (<fi 2 (x)) from the finite terms of the DeWitt-Schwinger expansion, and 
some concrete examples are given for scalar fields in four- and five-dimensional black hole spacetimes. The results are 
summarized in Sec.fV] In the Appendices lAl and [Bl we develop some formulas needed in the main part of the work. 

II. VACUUM POLARIZATION IN d-DIMENSIONAL STATIC SPACETIMES 
A. Green's Function Connection to ((f> 2 (x)} 

For a scalar field 4>(x) in a curved spacetime background we start with the action 

S = J d d x^/\g~\C, (1) 

where g is the determinant of the d-dimensional spacetime metric, and C is the Lagrangian for the scalar field </>, given 
by 

C=^{g^(x)4>(x),^(x), v - [m 2 + ^R(x)]^(x)}. (2) 

Here m is the mass of the field quanta, and it is assumed there is a coupling between the scalar and gravitational fields 
of the form £ R(x)4> 2 (x), where £ is the coupling constant and R(x) is the Ricci scalar of the background spacetime. 
Minimal coupling corresponds to £ = 0, while for £ = jjzj (£ — | m d — 4) the field is conformally coupled when 
m = 0, i.e. the action is invariant under conformal transformations of the type g^,v{x) —> g^vix) = ^(x) 2 g )iv (x) and 
<j)(x) — > (j)(x) — Q(x)( 2 ~ d ^ 2 tfi(x). Varying the action of Eq. ^ in relation to <fi gives the equation of motion for the 
field, 

(□ + m 2 + £R(x)) 4>{x) = . (3) 

This is a generalized covariant Klein-Gordon equation, where □ = g^ v V u — \g\^ x ^ 2 [(|g|) 1 ^ 2 g piy </ , i/] p is the 
Laplace-Beltrami operator in ci-dimensional curved spacetime. 

Quantization reveals that the field is composed of particles obeying certain commutation relations, and one wants 
to know how these particles move or propagate in the given curved background spacetime. A propagator is usu- 
ally defined by modifying the Klein-Gordon equation, Eq. ([3]), so as to include a source term J(x) such that 
(□ + m 2 + t;R(x)) <p(x) = J(x). This equation may be solved using the standard theory of Green's functions. One of 
the Green's functions, or propagators, that can be defined is G-p(x,y), which satisfies 

(a + m 2 +£R(x)) G F (x,x') = -\g{x)\' 1 ' 2 5 d {x - x 1 ) , (4) 

where 5(x) is the Dirac delta function. The solution for <f>{x) is then cj){x) — 4>o(x) — J d d x \g(x)\ 1 ^ 2 Gp(x, x') J(x'), 
where 4>o(x) is a function that satisfies the Klein-Gordon equation without a source term and (j>{x) corresponds to a 
quantum field operator acting on some state. 

Interestingly, vacuum expectation values of products of field operators can be identified with the Green's function 
of the wave equation, as we now show. The propagation of a free test particle in a vacuum |0) can be described by 
the correlation function G + (x 1 x') = (0\<p(x)(/)(x')\0} = (4>(x)(f>(x')}, where 4>{x') creates a particle at t' , which in turn 
is annihilated by cf>(x) at t. This makes sense if t > t' . Analogously, the correlation function G~(x,x') = (cj)(x')<f)(x)}, 
describes the propagation of a particle created by <f>(x) at time t, which in turn is annihilated by 4>(x') at t' . This makes 
sense if t' > t. To obtain a correlation function, or propagator, that has physical meaning in relativistic quantum field 
theory, either G + (x,x') or G~(x,x') is used, depending on the sign of the relative time. So to obtain a physically 
meaningful propagator that combines both G + (x,x') and G~(x,x') we can use (T (4>{x)4>{x'))) , where Dyson's time 



5 



ordering product operator T is defined as T (</>(x)</>(x')) = 6(t — t')<p(x)<p(x') + 9 it' — t)<p(x')(/)(x), with 9(t) = 1 for 
t > and = for t < 0. To call T (4>(x)4>(x')) the "time ordered product" is apt since the operators occurring 
under the symbol T are arranged from right to left with increasing times. Such a propagator is called the Feynman 
propagator, and one can show that this time ordered product of fields is indeed the Feynman Green's function defined 
by Eq. Q, i.e. 

iG F (x,x') = {T(<f>(x)ct>(x'))). (5) 

Using the Klein-Gordon equation, Eq. ([3]), and the properties of the step function 9 it' — t), one finds 

(□ x + m 2 + £R(x)) (T (4>(x)4>(x'))) = -i\g(x)\-V 2 6 d (x - x'). (6) 

Care should be taken since the step functions are time-dependent, and instead of zero the result is a distribution 
8 d {x — x') concentrated at equal times. Thus it follows that the vacuum expectation value (T i<pix)<j){x'))) is essentially 
one of the Green's functions of the covariant generalized Klein-Gordon operator, and we are justified in calling it the 
Feynman Green's function G-p(x,x'). In other words, the analysis shows that the Feynman propagator is a Green's 
function of the Klein-Gordon equation. 

Usually in quantum field theory the equation connecting Green's functions and expectation values, such as Eq. 
gives a bridge between the theory of propagators, in which scattering amplitudes are written in term of Green's 
functions, and the theory of operators, where everything is written in terms of the quantum field 4>{x). One finds the 
operators and expectation values, thus obtaining the Green's functions important for interaction theory. We see that 
in our study, the connection is inverted - we want (j> 2 ix,x') at the point x by expressing operator theory in terms of 
the Green's function and so we calculate the Green's function. Thus, Eq. ([5]) operates as a kind of duality. 

Since we are interested in the coincidence limit, Feynman's Green function is the best to use because it is more 
physical and also because the boundary conditions allow a Wick rotation of the equation to Euclidean space, where 

Gp(t, x; t', x') — —iGv{iT, x; ir , x ). (7) 

The Euclidean Green's function, Ge, now obeys 

(□ E - to 2 - £R(x)) G E (x, x ) = ~\gix)\- 1/2 S d (x - x') . (8) 

where De is now the Laplace-Beltrami operator in d-dimensional curved Euclidean space. There are advantages 
to working in Euclidean space. For instance, elliptic operators are more easily handled than hyperbolic operators, 
and after obtaining the Euclidean results one can Wick rotate back to Lorentzian spacetime using Eq. ([8]) since the 
boundary conditions for the Feynman propagator are automatically imposed by this procedure. 

The Feynman Green's function, or alternatively the Euclidean Green's function, is defined in terms of expectation 
values of products of field operators in the pure vacuum state. This is fine for describing the system at zero tempera- 
ture. To go further and describe a system at nonzero temperature one has to take into account that the system is no 
longer in a pure state, it is statistically distributed over all possible states. The full weight of statistical physics must 
be used, and the Green's functions are given by the average, suitably weighted, over all pure states of the expectation 
value of the products of field operators in those pure states (see e.g. Ref. (62j). 



B. Calculating the Green's Function 

The standard approach now used for calculating ((f) 2 (x)) was laid down by Anderson (40l |, based on earlier works 
by Candelas and Howard [2l|, [22j • We start with the Euclidean metric for a static spacetime in d dimensions with 
line element 

ds 2 = f(r)dr 2 + hir)dr 2 + r 2 dVt 2 . (9) 

Here r is the Euclidean time, r = —it, r is a kind of radial coordinate, and f2 represents a {d — 2)-dimensional angular 
space. The only restriction for this method is that the metric must be diagonal. The expectation value (<p 2 (x)) is 
found from the coincidence limit of the Euclidean Green's function Ge(x,x) = lmi x i^ x G~e(x, x'). For a scalar field 
in a spacetime given by Eq. ©, Ge(x, x') satisfies [see Eq. ©] 

(□ E - m 2 - £R(x)) G E (x, x') = ^=<S(r - r')<5(r - r')<5(tt, to'), (10) 
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where De is the Laplace-Beltrami operator of the Euclidean metric corresponding to Eq. ([9]) . Assuming a separation of 
variables, the independent homogeneous equations for r and f2 may be solved. Standard Green's function techniques 
then tell us that the r and f2 dependence of G-e(x,x') is equivalent to a representation of the corresponding delta 
function. We therefore use 



5(n,n>) = E E Y L{H} (n)Y* {H} (n>) , (11) 



as the Ansatz for the angular dependence of the Euclidean Green's function. The function Yi^.y in Eq. (|lip has been 
generalized to the set of hyperspherical harmonics. In four dimensions these are the usual spherical harmonics such 
that Ef E m ^.m (^)^e*m (^') = 3^X^(2^ + 1)-FH^ ■ ^')- As the Ansatz for the time dependence of the Euclidean 
Green's function, an integral or a sum representation is used depending on whether the scalar field is at zero or 
nonzero temperature, respectively. If the scalar field is at zero temperature, then 

5{t -t') = — / duoe luj{T - T > . (12) 

If the scalar field is at nonzero temperature T, then the Green's function is periodic in r — r' with period T^ 1 62], 
and a suitable representation for the delta function is 

oo 

&{t-t')=T exp [m27rT(r - t')] • (13) 

n— — oo 

Henceforth denote e = r — r'. If the scalar field is at zero temperature, then 

1 f°° 

G B (x,x') = — ckje^J2J2 Y ^M n ) Y U ^%X^(r,r'). (14) 

where x^(r, r ') IS the last component of the variable separated Green's function - a radial mode function. On the 
other hand, for a scalar field at some nonzero temperature T 



G E (x,x') = - J2 eIK£n EE^te}(")^}^')x^(r-,r'), (15) 

n=-oo i {n,} 

where k = 2-kT. In both cases the radial function obeys a differential equation obtained by putting the above expres- 
sions into Eq. (|10[) . Using this expression for the Green's function, and the preceding discussion on the connection 
between the Green's function and the expectation value (4> 2 (x)) allows us to calculate (<p 2 (x)) in the Hartle-Hawking 
vacuum. 

There are three difficulties when evaluating these expressions in the coincidence limit. The first is that the equation 
of motion for the radial function \ni (equivalently Xut) is quite complicated, with exact solutions only available for 
zero frequency. Asymptotic solutions are obtainable in closed form for massless fields on the horizon of Schwarzschild 
and Reissner-Nordstrom black holes [U, |4(| H3] • Partially analytical and numerical evaluations of the radial modes 
occupy the bulk of current research on this topic and will not be discussed here. The other two difficulties are that 
the sums over both i and n (equivalently the integral over to) produce divergences. The divergence resulting from 
the sum over t is actually only an apparent divergence and may be easily remedied. The standard trick is to realize 
that, given the delta function dependence, we are free to add a term proportional to the delta function. The large 
I contributions are then eliminated with the help of a WKB approximation. It is the divergence resulting from the 
sum over n which is more serious, and it is to this matter that we direct our attention. 



III. DEWITT-SCHWINGER RENORMALIZATION IN d DIMENSIONS 



A. General treatment 



To assign a physical meaning to ((f) 2 (x)) ^ it must be rendered finite via some renormalization process. The divergence 
resulting from the sum over n in Eqs. (I14|) and (|15p is related to the high frequency behavior of the scalar field. The 
high frequency modes of the field probe the spacetime geometry in a small neighborhood of an event. Since the metric 
changes negligibly in this neighborhood an adiabatic, short distance approximation for the propagator should give 
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the same divergent behavior as Eqs. (fT4| and (fT5|) . Isolating the ultraviolet divergences with such an approximation, 
these divergent terms can then be subtracted from Eqs. (fT4| and (| 1 5 [1 : leaving the renormalized, finite part of the 
Green's function. The now standard approach is to renormalize the expression for Gk(x, x') via the point splitting 
method of Christensen applied to the DeWitt-Schwinger expansion of the propagator (q.llol Il2l. [l3j. In d dimensions, 
the adiabatic DeWitt-Schwinger expansion of the Euclidean propagator is [13, Eq. 3.10] 



G» s (x,x') = 



7TAV2 



( 47 ™) d/2 ^0 



OO 



a,k(x,x') 



d 



dm 2 



2im' 



l-d/2 



H 



(2) 
d/2- 



(16) 



This equation is slightly different than that found in Ref . [131 ] , where the expression is given for the Feynman Green's 
function rather than the Euclidean Green's function and uses a different sign convention; the two are related by 
Eq. ([7]). Equation (fTB]) introduces several new variables that must be defined. Let s(x,x') be the geodesic distance 



between x and x', then define 2a(x, x') = s 2 (x,x') and z 2 = —2m 2 a(x,x'). The coefficients a,k(x,x') are generally 

(21 

referred to as DeWitt coefficients. The function H v (z) is a Hankel function of the second kind. Lastly, A(x,x ) — 
y/g(x)D(x, x')y/g{x') is the Van Vleck-Morette determinant, where g(x) — det(<7 M „(x)) and D(x,x') = dct(— cr-^). 
Expressing —2m and dm 2 in terms of z and dz (for fixed a), Eq. (|16|) can be written as 



Gg S (x, x') 



rA x /2 



( 47r )d/2 



^2a k (x,x')(-2 



m 2\(d/2-l-fc) z -2{d/2- 



1-fe) 



fe=0 



z<9z 



y d/2-lo-(2) 



By the derivative formula for Bessel functions [6S 







and defining ^ = d/2 — 1 — k, Eq. p?]) becomes 

-rDS' - -ittA 1 / 2 



G*e (a;, a;') 



(47r) d / 2 



^a fe (x,a;')(-2m 2 )^^i/( 2 )0 



(17) 



(18) 



(19) 



k=0 



The idea here is that the DeWitt-Schwinger expansion results from a WKB expansion for the Euclidean (or Feynman) 
propagator for a generic spacetime when the point separation is small. For a particular spacetime this procedure does 
not give the correct results for the Green's function with finite point separation, but it should reproduce the same 
divergent terms in the coincidence limit. Therefore, if the divergent terms of the DeWitt-Schwinger expansion can be 
isolated, then these will be the terms to subtract from G^(x, x') in order to make it finite 

(2) 

The Hankel function is related to the usual Bessel functions by Hy (z) = J y — %Y u {z). Note that z = i\z\ is purely 
imaginary in Euclidean space. For a purely imaginary argument one finds [63j 



H^(i\z\) = J v (i\z\) - iY u (i\z\) = i v I v (\z\) - i 
leading to 



i^i v (\ z \)-~(-iy Kv {\z\ 

7T 



2% v I v {\*\) + -i{-i) v K v (\z\), (20) 

7T 



2\v\ 



(47r) d / 2 



fc=0 



(-iy*i v {\z\)+%K v <\z\) 



(21) 



Since we are working in Euclideanized space the physical renormalization terms come from the real part of this 
expression, which will leave only the A"„(|z|) terms. The asymptotic behavior of .K„(|z|) for small argument, z —> 0, 
is 



KJ\z 




v < 0, 



In 



(22) 



clearly only those terms of the sum for which v > produce divergent terms in the coincidence limit. Let kd be the 
largest integer that is less than or equal to d/2 — 1; for even dimensions kd — (d — 2)/2 while for odd dimensions 
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kd = {d— 3)/2. It is clear that the divergences arise from those terms for which k < kd- Note however, that the terms 
for k > d/2—l are also divergent in the limit of vanishing field mass. For a massless field, where the DeWitt-Schwinger 
formalism is obviously not well defined, the mass must be replaced in those terms which diverge logarithmically as 
772 — > with a parameter \x which must be fixed by a renormalization condition, or perhaps, as claimed in (To| . must 
ultimately be determined experimentally. Taking the real part of , the divergent terms are 

G **M = ^^^a k {x,x'){2m 2 r\z\- v K v {\z\). (23) 



k=0 



The Van Vleck-Morette determinant A 1 / 2 and the DcWitt coefficients a k {x,x') both depend on the point separation 
z, or equivalently a(x, x') = ^a p cr p7 where a p = a ;p [llj. Essentially, a p is a vector that points from x to x' and has 
length equal to the distance from x to x' . Consequently, a p — > in the coincidence limit. 
Any scalar function may be expanded in a covariant Taylor series of the form [l4j 

= E ~TT V ^ ■ ■ ■ V<* k m<r ai ■■■° ak - (24) 



Christensen has calculated these expansions for A 1 / 2 , a (x,x'), ai(x,x'), and a 2 (x,x'), which are provided here for 
completeness [ill, Ell ■ He finds, with <Zq = 1, 

A 1 / 2 = i + ^R ap <r a ^ - ±R a ^a a a (T-i + {^gRafiRys + m RP Jp R fnr6 + ^R^s)^ ^ ^ a 5 + • • • , (25) 
a l = {\-(i)R-\{\-OR a a a 

gQRapR P ' p + TEo RP Rpar/3 + i^RprnaR P fj + l20 R af3;p P (40 _ 6^ R ' a ° 17 ~^ ' ' ' ' 



0.2 



ih RPTR p r + -m Ral3pT R ap P r + §(§ - C)R-p p + \{\ - 2 R 2 h — • ( 27 ) 



For calculations of (<j) 2 } up to d ~ 7, only the DeWitt coefficients up to a 2 are required for renormalization. For 
d = 8,9, a 3 is needed and has been found by Gilkey [16(. For d — 10, 11 a 4 is needed and has been calculated in 
the coincidence limit [13, El, EH . For higher dimensional spacetimes, subsequent a n coefficients must be calculated. 
Note however that for calculations of (T MJ y) or other quantities involving derivatives of the field, more of the a n may 
be required for a given dimension. Since A 1 / 2 and the a k contain powers of a p , the entire expression for Gdw should 
be expanded in powers of a p before taking the coincidence limit x — > x' . Let A 1 / 2 and a k be expressed as 

A 1 / 2 = aJ /2 + A; /2 + A 2 /2 + ---, (28) 

and 

a k = a° k + 4 + a\ + ■ ■ ■ , (29) 
where the j^ 1 numerical index indicates the corresponding term of Eqs. (|25p ~ (|27l) containing j powers of a p (note 

1 /2 

that A/ = 0). The bracket notation is the usual notation found in the literature indicating the coincidence limit, 
e.g. [afc] = a? k is the term containing zero powers of a p . 

These expressions can be put together and expanded to the appropriate order for any dimension. The result, 
however, would still not be in a form that can be combined with Eqs. (fTl| - (|15p . It would therefore be useful to 
workers in the field to have a compact formula for the renormalization terms as applied to calculations of ((j> 2 (x)) and 

The modified Bessel function .K"„(|z|) behaves differently for even and odd dimensions, so they must be considered 
separately. In even dimensions v is an integer while for odd dimensions v is a half integer. In the small z limit one 
may verify that Iz^if^dzl) behaves as 



= L 2 ^W(u- - ■ + ^TTTvTh^T ( > . ~ 2 ( * V + 7 ) ) + 0(\z\) (30) 

n = l ^ 
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for integer z/, and as 



i*r^(W) = E w 



^±2 (_ 1 )«+n+i 2 2n - 1 r(n - 1/2) (-1) 



I>-ri + 3/2)|z| 5 



2 l '+ 1 i> + 1) 



0(1*1) 



(31) 



for half-integral v. Incidentally, from these expansions for i£^(|z|) one can begin to see the connection with the 
Hadamard form of the Green's function (54|,[55|. These expressions imply that multiplying A 1 / 2 , a*, and "2f„(|z|) 
requires expansions of A 1 / 2 and a k to order 2u in a p prior to taking the coincidence limit. Some authors refer to this 
as the "adiabatic order." Using the expansions of a k and A 1 / 2 , Eqs. (|2"5|) and (|2T)|) , we may collect a^A 1 / 2 in powers 

of (jP, 



oo p 

a k A^ = a°A^ 2 + (alA 1 / 2 + a^ 2 ) + (a°A^ + a\A\' 2 + a 2 A^) + ■■■ = [a k ][A^] + £ £ a£ 

p=l j=0 



(32) 



The summand of Eq. (|23p may be expanded explicitly in powers of a p , giving 

2n p 



[a k ][A^](2my\zrK u {\z\) + (2m 2 f £ — ^ 1} " ~' T "" 



+ (2m 2 r^- 



+ 1 I>-n+l)|z| 2 " fc "-•> 

v ;l 1 p=l j=o 

(-i)"-"r(«) f> v a . A V2 

, 2^- 2 "+ 1 r(i/-n+l)U| 2 ™ ^ ^ fc 

1 v p=2n+lj=0 



(-1)" 



2"+ 1 I> + 1) 



2(lnY+7 



- 1 

E 1 



oo p 



p-i 



(33) 



j P =xi=o 



for integral v (even dimensions), and 



v+4 



[a k ][A 1 / 2 ](2m 2 r\zrK u (\z\) + (2m 2 )" E Wi 



(_ 1 )^+n+i 2 2n-l r ^ n _ 1 / 2 ) 



■J 2 I/ + 1 r(^-n + 3/2)|^| 



2n p 



a J A 1/2 

J;' 1 /< ./ 

p=l j=0 



v+i 



+ (2- 2 )^ E St 



(-l)»+ n +i2 2n - 1 T(n- 1/2) 



^ 2" +1 r(i/-n + 3/2)1 

n— 1 v 



z |2n-l k P-3 



p=2n+l j=0 



(-1) + 2?r /2 

2 v+1 T(v + 1) fc P-J 

\ ' p=l i=0 



(34) 



for half-integral v (odd dimensions). In the second term ctfcA 1 / 2 has been expanded to order 2v. Since a? k A x J 2 ^ is 
proportional to (a a ) p , it is clear that the third and subsequent terms all vanish in the coincidence limit, leaving 



Gdiv{ x , x ) 



E 

fe=0 



N(2mYkr%(N)+EEE 



" 2 ™ P 2 2 "- 1 (-m 2 )-"r(n) ajA^) 



n—1 p—1 j—0 



r(u-n+l) (pPo p ) n 



(35) 



and 

Gdiv^j a; ) = 



( 47r )d/: 



E 



fc=0 



2n P n2n-2(_ rn 2\u+n+\ V ( rl _ l\ J A 1 / 2 

[a k] [A^]{2m 2 r\zrK v {\z\) + e E E — • J ; 

n—1 p—1 j—0 



(<^p) r 



(36) 



r(i/-n + |) 

for even and odd dimensions, respectively; and where in the second term we have used \z\ 2 = m 2 a p a p . To reiterate, 
kd = (d — 2)/2 for d even and kd = (d— 3)/2 for d odd. While at this stage the even- and odd-dimensional equations 
appear to have the same form (with the simple replacement n — > n — i), it is clear from Eqs. (|30[) - (|31[) that the end 
result is not the same. In particular, as is known, the even-dimensional result contains a logarithmic divergence while 
the odd dimensional result does not. 

The preceding equations are covariant expressions that isolate the divergences in a generic d-dimensional spacetime. 
To perform any meaningful subtraction of these divergences from the Green's function, these terms must be expressed 
in a form commensurate with Eqs. (|14[) - (|15p . In particular, it would be nice if these terms could be expressed either 
as an integral over w or as a sum over n. It will be shown that useful integral and sum representations compatible 
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v Extra Terms 



_ / — c 4 

" 2 y 2c 2 ISJ 

24e 4 , 4 



_ (^ + i(2 C i- C 2 C6 )-^( 

2t[S|-^(^ + l|(^- 4C ^))+^ 



e * c 2 
3 

2 " 

192e 4 24 



^ + 4(3c!-2c 2C6 ; 



8 



10cJ - 12c 2 c 4 c 6 + 3clc 8 ) 







2c| 

2 2 



+ TSf-O* - 4c2C6 ) + I? ( 33c 4 " 36c 2 c 4 c 6 + 8c!c 8 ))] 



2 I rn c i 12 



(2c 2 , - C2C6 )) + | (5c! 



5c 2 c 4 c 6 + cic g 



)] 



~-~ ■ JJ T + if (3c! - 2c 2 c 6 )) + f| (lOcl - 12c 2 c 4 c 6 + 3c 2 . c 8 )] + if (35c| - 60c 2 c!c 6 + 20c|c 4 c 8 + lOclcf, - 4cici ) 



TABLE I: This table shows the extra terms generated by the modified Bessel function iCG-zD when one makes the replacement 
|z| 2 — > 2m 2 C2 n e 2n - Integer values of v are applicable to even-dimensional spacetimes, whereas half-integer values of v 

are applicable to odd-dimensional spacetimes. 



with Eqs. (| 14[) - (| 1 5 p can be found for even-dimensional spacetimes. At this time, however, a correspondingly suitable 
expression for use with odd-dimensional spacetimes remains elusive. Consequently, in what follows we primarily 
address renormalization with respect to even-dimensional spacetimes. 

Unfortunately it does not seem to be possible to obtain a simple, compact, general expression as an integral over u 
or sum over n. The first problem is that, while the second term of Eqs. (|35[) - (|36[) may simply be finite, as is the case 
for four dimensions, this is not generally true for higher dimensions, as will be shown explicitly for the six dimensional 
case below. These additional divergent terms may be addressed by Howard's method |25j, described in Appendix lAl 
and used below. 

As for the first term, we may proceed a little further but must use some care. Recall that the physical parameter 
approaching zero is e = t—t', then z must be expanded in powers of e with the end result that z 2 = —2m 2 Xm=i C2n,£ 2 ™ 
for some r-dependent coefficients cm- The C2n are combinations of the metric functions / and h, and their derivatives. 
Expanding z~ 2n in powers of e one gets a series of terms proportional to e~ n , e~ n+1 , . . .e _1 plus a constant term. 
This means that 

(2m 2 ) v \z\- v K v (\z\) = } 2m )l K v {mEy/J) + Extra Terms. (37) 
(rneV/r 

The extra terms, which will be denoted E v , must be determined for each and so far a compact expression giving 
the extra terms for a given v is unavailable, but may possibly be found from a lengthy exercise in combinatorial 
gymnastics. The extra terms for the first few integral and half-integral v are presented in Tables [J and expressed in 
terms of the coefficients C2«. 

In practice the extra terms are straightforward to calculate using a computer algebra system. One simply takes the 
difference of (2m 2 y\z\- ,J K^(\z\), with \z\ 2 -> 2m 2 J2n=i c ^e 2n and expanded around e = 0, and (2m 2 y\z\~ v K v (\z\) 
with the replacement \z\ 2 — > 2m 2 C2S 2 . One important point that should be noted is that finding the coefficient C2 n 
in the expansion of \z\ requires one to first calculate to order e 2 ™ -1 . In Table HT1 the extra terms are presented 
explicitly in terms of the metric functions / and h for the first few integral and half-integral values of v. In four 
dimensions k c i = 1, in which case v ranges from to 1 and these extra terms contribute no new divergences. In six 
dimensions — 2, v ranges from to 2, so one extra divergent term arises when v = 2. In eight dimensions the extra 
divergences come from both v = 2 and v = 3. Obviously a similar situation occurs for odd dimensions. 

All these divergent terms lurking about within Eqs. (|35[) - (|36[) must now be expressed as an integral over u or a 
sum over n, commensurate with Eqs. (|14p - p5D . To this end we use an integral representation of the modified Bessel 
function, some identities proved by Howard J2a| , and the Plana sum formula [Z(J EH, EE Efi| , 

f(n)dn = /(») - - * / ^F-T W + ^ M - l ^ > ( 38 ) 

Jj n=j L Jo e r 

to convert between integrals and sums. 
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Extra Terms (E„) 






f' 2 

2iph 



16e/ 5 / 2 h 



f /2 r 2 
V^f' 2 ( 

Sf s / 2 h ^ 



£ 3 / 



I 1 f 23/' 2 
60h ^"iw/ 2 " 



+ 



f 2 h 



+ 



2fh f )\ 

1 i 1 ) i ™ 2 *' ( A. I 7 f i 



4fh 



17/' , hi 
16/ T 2h 



£ f 23/' 2 , /V , h' 2 \ , £_ > h!_ _ f'\\ 

12 ^ 360/ 2 ~ 5/fi T 16h7 y ~ 6 ^ 5/ 8/1 8f J ) 



12 V 



23/' , V. _ 
24/ 1" 2fe /' 7 ) 



TABLE II: This table shows the extra terms generated by the modified Bessel function if„(|2|) for the specific metric of Eq. 
©■ 



B. DeWitt-Schwinger Renormalization in Even dimensions 

1. Renormalization formulas at zero and nonzero temperatures 

For v an integer it is shown in Appendix [Bl that an integral representation of K v {z) for small z is 

*"W = TTT^ilY I™ dtcos{zt){t> + iy-V\ (39) 

r(^ + 5 ) ^2/ j 

For T = one has to connect Eq. (f3"9"|) with the e l " e dependence of Eq. fT4")) . Consider the change of variables 
t = uj/y/rn 2 /, and z = mey/f; then 

-^^J^rneV?) = ; K^-^ /" ™{e^ 2 + m 2 /)- 1 ^. (40) 
{mey/f) v (-/) r( I/ + 5) 

This result generalizes the integral representation found by Anderson (23, Eq. (3.4a) and (3.4b)]. 

For nonzero temperature T, one has to connect Eq. (f3"9")l with the e 4Ke ™ dependence of Eq. (|15[) . We instead make 
the change of variables i = xk/ yj m 2 f to first obtain 

^H^-K^msy/f) = k^-u r cos( K£ *)( K V + m 2 /)"" 1 /^. (41) 

The Plana sum formula, Eq. (|38[) . enables the integral in this equation to be converted into a sum plus some 
residues and is valid if the function / satisfies three conditions: (i) /(r + it) is holomorphic for r > j for any t, 
(ii) lim^oo \ f(r + it)\ e~ 27r l f l = uniformly for every r > j, and (hi) lim T ^oo dt |/(r + e 27r ' * I = 0. A naive 
application of the Plana sum formula would be to use j = 0, corresponding to the lower limit of integration in Eq. 
(|41j) . However, for ,7 = the integrand of Eq. (l4Tj) is not holomorphic at r = 0. Consequently, one must break up the 
integral into two parts 

dx cos{nex){n 2 x 2 + m 2 fy- l/2 = / dx cos{nex){n 2 x 2 + m 2 ff- l/2 + / dx cos{nex){n 2 x 2 + m 2 fY~ l/2 . (42) 



For the first integral cos(nex) w 1 in the coincidence limit and the solution may be expressed as a hypergeometric 
function depending on v [63, Eq. (2.271)], 

jf 1 dxcos{ K ex){n 2 x 2 + m 2 /)^ 2 = (m 2 /)"^^ Q, \ - u, | . (43) 

In general, this hypergeometric function is equivalent to a polynomial in half integer powers of (k 2 + m 2 /) plus a 
logarithmic term. Applying the Plana sum formula to the second integral gives 

cos(k£x)(k 2 ie 2 + m 2 f) u ~ 1/2 dx = ^ cos(Ken) (k 2 u 2 + m 2 /)" 2 - -(k 2 + m 2 f) v ~^ 

n=l 



r 
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Putting this together, we have 
(2m 2 Y 



(mevT)" 



K v {me^/J) 



r- < cos(Ken) (k 2 



n 2 + m 2 /) ~ 2 -^ 2 + ™ 2 fY^ + (m 2 /)"-^ r, o " "> o> 



11 3 k 2 



2' 2 '2' m 2 / 



{ [(1 + 4<) v + m2/] _ [(1 _ zt) v + m2/] - 2 }} (45) 

This result generalizes the sum representation found by Anderson [40j, Eq. (3.7a) and (3.7b)]. Finally, the renormal- 
ization terms for the d-dimensional spacetime of Eq. ([9]) are 



(4^/2 Z- 



E 



(-/)T(i/ + |)7 



3(we)(w 2 +m a /) , '- 1 / a dw 

!/ 2n p 



n— 1 p— 1 j—0 



2 2 "~ 1 (-TO 2 )^"r(n) ojAp^- 
r(z/-n + l) (cr^ffp)" 



(46) 



for the case of a scalar field at zero temperature T = 0, and 



Gdiv^j a; ) 



E 



(^^1 (_/)^r(i/ + |) 



k=0 

— i 



cos(Ken) (n 2 n 2 + m 2 f) v 2 - -{k 2 + m 2 /)" - * 



{ [(1 + it) V + m 2 /] ^ 1/2 - [(1 - it) V + m 2 /] ^ 1/2 } 



11 3 k 2 



t 2' 2 "'2' m 2 / 
for a scalar field at nonzero temperature T > 0. 



i/ 2n p 2n-l/ ~,2\i/— n 



3AV3 



^ + EEE- 

n— 1 p— 1 j— 



(-m 2 )-"r(n) 

]?(z/-n + l) (crP(Tp) r 



(47) 



2. Examples: d = 4 and d = 6 



The formulas given above provide simple expressions to calculate the renormalization terms for the generic even- 
dimensional spacetime of Eq. @. Below we mention the case d = 4 and study more carefully the case d = 6. For any 
d-dimensional spacetime with line element given by Eq. ([9]), we generalize u M to fill. |42| 



'\2 



24 fh 



e 5 ( f' A 3 (f'yh' 3 / /2 /" \ r 
120 I 8f 2 h 2 + 16 //i 3 8 fh 2 ) + (£ ' 



(48a) 



Ah 



5" 

24 



8h 3 Ah 2 



(48b) 



^=0 i = l...d-2. 
and (|4"7) . we use the values for E v as listed in Table ITT1 



(48c) 



In applying Eqs. 
d = 4: 

It is straightforward to show that in four dimensions Eqs. (|4"6")l and (|4"T|) are identical to those obtained by Anderson 
pol. l42l] and used by several subsequent authors. After letting e^Owe find, for T = 0, 



G div (x, z') = -^LJ™du; j(u; 2 + m 2 ff' 2 + ~ ^ - J2 



'V + m 2 /r 1/2 



192/ft V r 



2T 

/' 



2/' 

T 



77 



(49) 
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One may verify that expanding Eq. (|47|) correctly reproduces the results obtained by Anderson for T > 0. 
d = 6: 

Another, less trivial, example can be given for a scalar field in six dimensions. We consider three classes of spacetime: 
spherical, flat, and hyperbolic, corresponding to & = 1, 0, or — 1, respectively. Consider first the last terms of Eqs. 
(|46f and (|47|) . Using the values of E v given in Table HT1 and calculating the last sums, we find 



E 

k=0 



v 2 ™ V 2 2n-l 



n— 1 p— 1 j—0 



r(i/ - n + 1) (Vo-p)" 



e 2 6f 2 h \r f f ft 1 



where 



(50) 



eft = - 2 



24/ft \r 



2/" 2f 



Ti 



JLil 

fh 2 [ r 3 

f 2h r 



5/' 2 /_9_ 
6/ 2 ^25 
ft" ( 1 



12 / 2ft V25 



/' uo 



3ft V 20 



3/ft V5 V ' 6ft 2 \50 "J ' f> \f \180 
f" \ f ff' 2 [27 \ Ifh! A27 A 5h' 2 

f \8f 2 \A0 V 48/ft ^210 V 48ft 2 V 100 



/" (If A 27 



/ V24/ V210 



15/' 



03 



5/i' / 61 



24ft V 100 
7ft' 3 



r 



ii 



12/' \20 
19/"ft' 2 13ft'ft" 



/' A3 



12/ V20 ^ 
/"ft" ft'/'" 



2ft ~7~ 

ft'" /< 4 ) 



240ft 3 480/'ft 2 480ft 2 60/'ft 40/'ft 240ft 120/' 
Using Eq. (|A5j) for the integral representation of e -2 , the T — divergent terms in the limit e — > are 



(51) 



Gdiv(^, X ) 



(47T) 



do; 



^( W 2 + m 2 /) 3 / 2 + 2(£ - I)i?(c 2 + m 2 ff' 2 + [a 2 ](^ 2 + m 2 ^ 2 

f'u (8 2/" 3/' ft' 
6/ 2 ft /' / ft 



(52) 



where, for d = 6. 



24^tS 7i(#ft 



1) [(-41 + 420£ - 1080£ 2 ) + £ft (29 - 360£ + 1080£ 2 )] 

+ w {fr ~ t) K 29 " 29 °£ + 720 £ 2 ) " 5kh ( 5 ~ 54 £ + 144 £ 2 )] 

108(1 - 60 2Z ^ + T " 2 r") + ( 23 _ 140 £ + 120 ^ 2 ) /=" ~ 4 ( 



ft'- 
TP 



4 - 85£ + 330$ 

') (-15 + 56£ + 96£ 2 ) ^ + 4 (1 - 40$ + 180£ 2 ) ^ - 40(-l + 5£)^ 
+ i f- (23 - 140£ + 120£ 2 ) ^ + (-33 + 140£ + 120£ : 
+ 2^ ((21 - 130£ + 120£ 2 ) £j - (-13 + 40£ + 120£ 2 ) £ 

- 4 d - m (¥ (ef - 1 4 ^) - ^ ( 4 £ - 13^) + 4 ^ - 2Z ^) 

5 (5 - 26£ + 6£ 



2 \ 



21 - 110£ + 30£ 2 ) ^ 



70£-30£ 2 )^ 



off 

2(1 -r,;,(l4^ 



( 2 < 



13 - 70£ + 30^ 2 ) £ 



2 (-13 

5 (5 - 27^ + 12^ 



2\ /i'- 



h X* f 



13 



+ 8^ — 



2\ W_ 

h 



10 (1 - 6f + 6^ 2 ) C 



7~ ( 5 7 



:h' 



z / ft 



(53) 
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On the other hand, for nonzero temperature T > 0, Eq. (|Aip may be employed to find a sum representation of e 2 
and one may show that 



Gdiv(x, x ) 



(An) 3 



E 



3^(«V + m 2 /)3 /2 + 2 ^ _ l )i?(K 2 n 2 + m 2 /)V2 + ^ ^2 + m 2 /} - 1/2 



/8 2/" 3/' 



/' / 



In 



« + \/ k 2 + m 2 f 



im 4 + (e-i)m 2 i?+M 
2 o 



4ift f 00 
37^ 7o e 



- { [(1 + ^) 2 ft 2 + m 2 /] 3/2 - [(1 - ^) 2 ft 2 + m 2 f] 3/2 } 



= 27Tt 



T { [(1 + *t)V + m 2 /] 1/2 - [(1 - it) V + m 2 f] 1/2 } 



-[a2]^° ;5 ^ T {[(l + i t) 2 ft 2 + m 2 / ]- 1/2 -[(l-zt) 2 ft 2 +m 2 /]- 1/2 } . (54) 



These results are quite general and unwieldy, but expressions such as Eq. (|53[) . for [02], may become quite simple 
for particular situations. For a minimally coupled scalar field in a six-dimensional, asymptotically anti-de Sitter 
Reissner-Nordstrom black hole spacetime with 



3 r 3 r 6 ' 



(55) 



[02] reduces to 



N = 



4(14Q 2 - 5Afr J ) 



3\2 



75r 



16 



(56) 



a remarkably simple result (note that in six dimensions a minimally coupled field has £ = 1/5). 

Put in perspective, we have presented a formula for the renormalization terms of (<f> 2 (x)}, which in the particular 
important case of a minimally coupled scalar field in a (d — 6)-dimensional asymptotically anti-de Sitter Reissner- 
Nordstrom black hole spacetime, yields an extremely simple compact formula. The dimension d = 6 is of importance 
because it is the simplest even higher dimension that can be made compatible with the extra large dimension or brane 
world scenarios. 



IV. ESTIMATE OF ((j) 2 (x)) FOR MASSIVE FIELDS 



Using the DeWitt-Schwinger expansion we have isolated the divergent terms in the coincidence limit. These were 
precisely the terms of the expansion up to k = kd- The DeWitt-Schwinger expansion does not give the correct results 
for (4> 2 (x)), even after removing the divergences, because the expansion depends only on the local structure of the 
spacetime whereas the true field modes also depend in part on the global structure of the spacetime, for example, 
the effective potential around a black hole. However, when the field is massive enough, the higher order terms in the 
expansion may be used as an approximation to the renormalized value of (4> 2 (x)) such that 

(^(x)) » lim G Ten (x,x') = km 7-4^ E &' 3 a k (2m*y\z\- v K v (\z\). (57) 
x'^x z^o (Air) > z — ' 

V ' k>k d 

For k > kd, v < 0, in which case \z \~ v K v (\z\) = 2- {lJ+ ^T(-v) + 0(z 1 ) for small z in both even and odd dimensions. 
It follows that 

- 7^72- E i^n-v) 

y ' k>k d 

in the coincidence limit. 



(2 d TT d / 2 m 2 )- 1 [a d } + --- deven, 
( 2 V d - 1 )/ 2 m)- 1 V^i]+--- rfodd, 
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FIG. 1: Plot of (0 2 (s)> in a four-dimensional Reissner-Nordstrom black hole spacetime with mM — 2. |(a)| Near the horizon. 
From top to bottom the curves correspond to the cases \Q\/M = 1.0, 0.99, 0.95, 0.8, 0.0. The radial coordinate s = r/M — 1 — 
(Q/M) 2 is zero at the horizon. 



(b) On the horizon as a function of \Q\/M. 



For d = 4 one finds (<j) 2 (x))d=4 ~ (47rm)~ 2 [a2\- This agrees exactly with the result reported by Anderson |40l . 
Eq. (4.1)], but note that Anderson's equation contains a typographical misprint in the first term, which should be 
g(C — \)R;p' p ■ This misprint has no consequence since R = in the spacetime considered by Anderson, but would be 
important elsewhere. Calculating the coefficient [02] for a Reissner-Nordstrom black hole spacetime, for which 



2M 
r 



Q 2 



gives 



02 



45r* 



■(13Q 4 - 2AMQ 2 r + 12MV), 



(59) 



(60) 



leading to the near horizon behavior of {4> { x ))d=i plotted in Fig. 1(a) It may be seen that this correctly reproduces 



Fig. 3 of Ref. [40j. Figure 1(b) shows the behavior of (</> 2 (x)) d=A on the horizon as a function of the charge-to-mass 
ratio, where it can be seen that the value of (<ft 2 (x))d=4 increases rapidly to a finite value as the black hole approaches 
extremality. 

In Refs. [13, IHJ it is emphasized that the finite terms of the DeWitt-Schwinger expansion give a good estimate 
for (<f> 2 (x)) rcn when mM > 1, especially near the horizon. Since the horizon radius obeys ~ M and the Compton 
wavelength associated to m is A ~ 1/m, the rough inequality can be translated into Th/X > 1. This could have been 
expected on physical grounds. On one hand, particles with much longer wavelengths (lower mass) are outside the 
validity of the approximation and it cannot give good results. On the other hand, since vacuum polarization happens 
most intensely near the horizon, then particles with wavelengths on the order of the horizon radius or less (77,/A > 1) 
are well described by the approximation because they fit within the most probable characteristic geometric length of 
the fully quantum processes in a neighborhood of the horizon. 

d= 5: 

In five dimensions the metric function f(r) for an asymptotically flat Reissner-Nordstrom black hole is 

2M Q 2 



1 



from which we find 

N = 



(61) 



(62) 



[48MV + Q 4 (-17 + 460£ + 60£ 2 ) + 24QV [M - 30M£ + 2r 2 (-l + 5£))] • 

Notice that, unlike in four dimensions, in five dimensions [02] depends on the coupling constant £. From Eq. (|6ip it 
is straightforward to locate the outer horizon at r\ = M fl + \J\ — (Q/M) 2 ^j . Letting s = r 2 — r 2 , the near horizon 
behavior of {<fi 2 (x))d=5 is plotted in Figs. 2(a) and 3(a) for £ = and £ = 3/16 (conformal coupling) respectively. 
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Plot of ((j) 2 (s)) for a minimally coupled field (£ = 0) in a five-dimensional Reissner-Nordstrom black hole spacetime with 
Near the horizon. From top to bottom at s = the curves correspond to the cases \Q\/M — 1.0, 0.99, 0.95, 0.8, 0.0. 

r 2 /M - 



FIG. 

mM 

The radial coordinate 



a 



y/l — (Q/M) 2 is zero at the horizon, (b) On the horizon as a function of \Q\/M. 




FIG. 3: Plot of (4> 2 (s)) for a conformally coupled field (£ = 3/16) in a five-dimensional Reissner-Nordstrom black hole 
spacetime with mM = 2. |(a)| Near the horizon. From top to bottom at s — the curves correspond to the cases 
\Q\/M = 1.0,0.0,0.8,0.99,0.95. The radial coordinate s = r 2 /M - 1 - y/l - {Q/M) 2 is zero at the horizon, [(b)] On the 
horizon as a function of \Q\/M. 



The behavior of the conformally coupled field is quite interesting and different from the behavior in four dimensions. 
The value of (4> 2 (x)) on the horizon is plotted in Figs. 2(b) and 3(b) for £ = and £ = 3/16, respectively. The main 
feature is that while (4> 2 (x))\ x=Xh , the value on the horizon, increases monotonically with \Q\/AI for the minimally 



coupled field as for the four-dimensional case of Fig. 
which occurs for a subextremal black hole with \Q\ 



1(b) 



M 



for the conformally coupled field there is a minimum value 
.97. These results show that {(f> 2 (x)} is well behaved near 



the horizon as might be expected, but some interesting variation arises for a non-minimally coupled field. 



V. CONCLUSIONS 



In this paper we have derived a compact expression for the DeWitt-Schwinger renormalization terms in d even 
dimensions. Beginning with the general d-dimensional formula for the DeWitt-Schwinger expansion, the divergent 
terms in the coincidence limit were isolated by considering an expansion in a. From the properties of modified Bessel 
functions a useful integral representation for K v (z) in even dimensions was found in the coincidence limit. This 
integral may be used for the case of a scalar field at zero temperature without further modification. For a scalar field 
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at nonzero temperature T the Plana sum formula was used to convert the integral into a sum plus residual terms. 
The resulting formulas, Eq. (|46|) for a scalar field at zero temperature, and Eq. (|47|) for a scalar field at nonzero 
temperature T, are given in a form compatible with calculations of (cj) 2 (x)} and (T M „(a:)) in static spacetimes. These 
formulas will be particularly useful for calculating (<fi 2 (x)} and (T M!/ (cc)) in arbitrary black hole spacetimes of even 
dimension. The formulas found reproduce directly in the case d — 4 the results obtained by Anderson [4(|. As a 
further example, the renormalization terms were calculated for six dimensions. Christensen remarked that calculating 
quantities such as (T U y(x)) in dimensions greater than four "would be extremely long and would probably have to be 
done on a computer" [131 ] mainly due to the complexity of the renormalization problem. While the derivation of Eqs. 
(|46|) and Eq. (|47|) did not require any special computing power, it is certainly true that calculating quantities such 
as [a,2\ and for a particular spacetime would be quite lengthy and extremely prone to error without the aid of a 
computer. Lastly, the finite terms of the DeWitt-Schwinger expansion that are nonvanishing in the coincidence limit 
may be used as an approximation to (4> 2 {x)). It is shown that this reduces to a sum over the DeWitt coefficients, and 
is discussed in more detail in four and five dimensions. 

As we have emphasized, in classical general relativity Einstein's equations relate the spacetime curvature to the 
distribution of classical matter as encoded in the stress-energy tensor. Unfortunately the Universe is not so simple, 
being composed of quantum, rather than classical, matter. While some argue this indicates the need for a complete 
quantum theory of gravity, a first step, used here, is semiclassical general relativity, where the stress tensor appearing 
in Einstein's equation is replaced by the expectation value of the stress tensor of quantum fields. Despite this objection, 
semiclassical general relativity has provided us with some impressive results and deep insight into the behavior of 
the Universe, but calculating the renormalized expectation values of {4> 2 {x)) and (T^(x)) is quite difficult in curved 
spacetimes. We have dealt with this difficulty here for (i-dimensional static spherical symmetric spacetimes in even 
dimensions. However, it is clear that, as it stands, the semiclassical theory is inadequate as a complete theory of 
gravity at the quantum level. By using the expectation value of the quantum fields, information about fluctuations 
of the fields, a defining characteristic of quantum field theory, is lost. At the very minimum, the semiclassical theory 
must be extended to incorporate some notion of fluctuations, and work in this area is being done by several authors; 
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APPENDIX A: GENERALIZATION OF HOWARD'S IDENTITIES 

For e~ 2 , Howard [25| proved that 



for small, nonzero e. Using Howard's procedure for p an even integer this result is easily generalized, for the problem 
at hand, to 




Acknowledgments 




(Al) 




(A2) 



where B p are the Bernoulli numbers. For p an odd integer, the identity is generalized to 




(A3) 



An integral representation was proved by Anderson, et. al. [42j by noting that 




as 



e 



0. 



(A4) 
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By repeatedly taking the derivative of both sides of Eq. (|A4[) and then letting A — > 0, one finds 

1 _ (-1)™ r°° 



e 2n T(2n) J Q 
for even powers of e, and 

1 (-1)" +1 r°° 



dtt Zn - L cos(et) (A5) 



1 roc 

— / dtt 2n sin(et) (A6) 
1) Jo 



for odd powers of e. In Appendix [B] we must evaluate 

dtt 2n ~ l cos{et). (A7) 



This is done by writing 

[°° dtt 271 ' 1 cosfet) = [°° dtt 2 "' 1 cos(et) - [ dtt 2 "- 1 cos(rf) = (~ 1 )" r ( 2 ") _ J_ 
Ji Jo Jo £ 2n 

where the l/2n comes from expanding the solution of the second integral near e = 0. 



(A8) 



APPENDIX B: DERIVATION OF EQ. ([39} 

To obtain formula ([39]) . begin with the recursion relation for the modified Bessel function, 

K v+1 {z) = -K v {z)-K' v {z), (Bl) 

z 

and the integral representation for Kq(z) [63j, Eq. (9.6.21)] 

= f<^- < B2 > 

From the recursion relation K\(z) = —K' a {z). Taking the derivative of Eq. (|B2j) and integrating once by parts gives 

tsm(zt)dt , , . . o si io 



Ki(z) = / ; 9 1 ' 9 = lim sin(zi)(i 2 + l) 1 / 2 - z / cos(zi)(* 2 + if^dt. (B3) 
Jo (* + 1) 1 s -*°° o Jo 

Taking the limit z — ► 0, the first term vanishes and only the second term remains. Applying the recursion relation 
again to find K^iz), a fortuitous cancellation and another integration by parts results in 

K 2 {z) = Z 1J~ cos{zt){t 2 + 1) 3/2 ^ « (|) 2 ^ cos(z0(i 2 + l) 2 - 1 / 2 ^. (B4 ) 

By induction one is lead to 

K ^) = TTT^Q 2 I™ ™{zt){t 2 + ir^dt. (B5) 

T(v+ ij V2/ J 

One might worry that it is unreasonable to assume that the first term on the right hand side of Eq. (|B3[) vanishes, 
as the limit z — > should only be taken at the end of the calculation, in which case the additional z^ 1 multiplying 
Ki(z) leaves us with a problematic lim s _ +00 z" 1 sin(zt)(t 2 — l) 1//2 |g. 

In fact, we can check that this integral representation reproduces the correct limiting behavior of z~ u K v {z) when 
z goes linearly to 0. First let 

dtcos(zt)(t 2 + 1)"- 1/2 = / dtcos(zt)(t 2 + 1)"-^ 2 + dtcos{zt){t 2 + IY- 1 ' 2 . (B6) 

n Jo Jl 
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For z in the coincidence limit we may expand the first integrand, cos(zt)(t 2 + 1)" 1 I 2 w (t 2 + 1)" 1 / 2 , and solve the 
first integral 

(B7) 



(B8) 



dtcos(zt)(t 2 + I)"- 1 ' 2 = 2 Fi (|, | - i/, |, -1) + / dtcos(zt){t 2 + If- 1 ' 2 . 
For the second integral, expanding the integrand for t > 1 gives 



/oc oo 
dtcos{zt){t 2 + \y- 1 / 2 =T(v+\)Y^ 
n=l 



^ r(n)r (i/-n+§) A 



<ii cos(z<)t 



2 1/— 2ra+l 



This expression may be analyzed with the use of some identities proved by Anderson, et. al. [42j |. based on identities 
proved by Howard [25| , and discussed further in Appendix El It follows from the results of Appendix [A] that 



dtcos(zt){t 2 + ly- 1 ' 2 = r (v + i) \J2 



i 



jr(n)r(!/-n + §) 



(_l)"-n+ir(2i/-2n + 2) 1 



^2y-2n+2 



1 J 

rTw? (lnz + 7) + ^r( n )r(,-„ + f j ,, 



2(// -n + 1) 
^cos(zi)t 2,y - 2 " +1 1 . 



r(j/ + i)v^ 

Renaming indices on the last sum, rearranging terms slightly and combining with the first part of Eq. (|B7I 

(-iy- n+1 T(2v-2n+2) 1 



(B9) 



/ dtcos(zt)(t 2 + iy-^ 2 = r(i/ + i) 

"'O U=i 



E 



r(^ + i)v^F 



E 



r(n)r (i/ - n + 1) z 2 "- 2 «+ 2 I> + 1)0F 

r(z/ + i)0F 



lnz + 7 



^^-n + W^-n + f) ^2nT(i/ + n + l)r(|-n) r (i/ + §) 



r( " +1) # 2 F 1 a,i-,,i,-D 



i2' 2 '2' 



(BIO) 



The last sum may also be expressed in terms of a hypergeometric function, and the limiting behavior is finally 

r(^ + i)0F 

{2(u-n+l)T(n)T (f-n+§) 



to = V (-1) 2 — +1 r(2^-2n + 2)^F _ (-1)" 
1/1 ' f^ 1 2 v T(n)T{v-n+f ) z 2l '- 2n + 2 2 v T{v+l) 



In z + 7 + 



rfr + i) 
4(r(^ + 2) 



i F 2 (l,l,f,2,2 + ,,-l)- r ^ + ; ) 1 f 2 F 1 (l,l 



r(, + i) 



0{z). (Bll) 



This result does not look the same as the expansion given in Eq. (|30[) . but we checked the expansion explicitly and 
believe the agreement is exact. To have exact agreement, the last three terms in brackets of Eq. (|B11[) must be 
identical to — 5Z n=1 n~ x — In 2. The two expressions were evaluated numerically with a precision of 100 digits up to 
v = 100 and in each case were found to agree within the working precision. While this numerical evaluation does not 
constitute a rigorous proof, it is a strong indication that the two expressions agree exactly and therefore the integral 
representation is valid. 
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